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ABSTRACT 

The beaming effect (aka Doppler boosting) induces a variation in the observed flux of a luminous 
object, following its observed radial velocity variation. We describe a photometric signal induced by 
the beaming effect during eclipse of binary systems, where the stellar components are late type Sun- 
like stars. The shape of this signal is sensitive to the angle between the eclipsed star's spin axis and the 
orbital angular momentum axis, thereby allowing its measurement. We show that during eclipse there 
are in fact two effects, superimposed on the known eclipse light curve. One effect is produced by the 
rotation of the eclipsed star, and is the photometric analog of the spectroscopic Rossiter-McLaughlin 
effect, thereby it contains information about the sky-projected spin-orbit angle. The other effect is 
produced by the varying weighted difference, during eclipse, between the beaming signals of the two 
stars. We give approximated analytic expressions for the amplitudes of the two effects, and present a 
numerical simulation where we show the light curves for the two effects for various orbital orientations, 
for a low mass ratio stellar eclipsing binary system. We show that although the overall signal is small, 
it can be detected in the primary eclipse when using Kepler Long Cadence data of bright systems 
accumulated over the mission lifetime. 
Subject headings: stars: binaries: eclipsing 



1. INTRODUCTION 

The angle between the spin axis of each star in a bi- 
nary system and the orbital angular momentum axis 
is an important clue for studying stellar binary forma- 
tion and the effect of tidal processes on their evolution. 
Naively one would expect the spins of the two stars and 
the orbital angular momentum to be well aligned, as 
they all originate from the angular momentum of the 
same primordial molecular cloud. However, this sim- 
ple view might be misleading. For example, binary 
stars with high o bliquity could result from elongated 
primordial clouds (Bo nnell et al.l I1992T) . or cha otic ac- 
cretion during star formation ([Bate et aTl I2010T ). This 
could also be the result of interaction with a third body 
throu gh the Mazeh-Shaham effect (jMazeh &; Shahaml 
1979), which is ess entially a com bination of the so-called 
Kozai o scillations dKozail H~ 962) followed by tidal dissi- 
pation (jFabrvckv fc Tre mainc 2007). A large sample of 
binary systems with measured obliquities is needed for 
studying the above processes. 

To date, spin-orbit alignment has been measured, with 
spectroscopy, in only three eclipsing stellar binaries, two 
of which are aligned (V1143 Cyg jAlbrecht et"aHl2007t 



NY Cep, Albrecht et al 



(DI Her, lAlbrecht et al 



I20TT 
l200" 



and the third misaligned 
It is interesting to note 
that for transiting planetary systems the spin-orbit angle 
is bei ng measured for an increasing number of systems 
(e.g., IQuelozet al1l2000t I Winn et all [2001 ITriaud et all 
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I2010t iWinn et ail 1201 ID . currently including one third 
of the ~100 known transiting systems. Almost half 
were identified as not well aligned, thereby challenging 
the simplistic understanding of the formation and evolu- 
tion of planetary systems (e.g.. iFabrvckv fc Winnl [20091: 
I Triaud et aTJl2010l : (Winn et al.ll20l"ot see also lSchlaufmanl 
2010). 

We show here that the beaming e ff ect (aka Doppler 
boosting: e.g.. iZucker et al.l 120071: iMazeh fc Faiglerl 
120101 : iShporer et al.ll2010f ). combined wi th the R o ssiter- 
McLaughlin effec t (herea fter RM; iRossiterl 119241 : 
iMcLaughlinl 119241 : see also lHoldli893 ). induces a pho- 
tometric signal during eclipse, where the shape of this 
signal is sensitive to the spin-orbit angle. Therefore, cor- 
rectly modeling the light curve allows to measure this an- 
gle. To t he best of our know ledge this approach was first 
used bv iHills fc Paid (|1974f ). for a wh ite dwarf. M ore 
recently, this approach was noted byjLocb (200 5)) for 
transiting planetary systems, and bv lvan Kcrkwiik et al.l 
(2010) for eclipsing w hite dw a rfs. T his signal was the- 
oretically modeled by iGrootl ([201 ll ) for three kinds of 
syste ms: eclipsing double white dwarfs (fStcinfad t et al.l 
l2010f K eclipsing O type binaries, and transiting planets 
orbiting an early type star. Here we present a further, 
fourth case: stellar binaries whose components are Sun- 
like stars, of spectral type F6 and later. Those stars are 
fundamentally different than earlier type stars and white 
dwarfs, as they are slow rotaters. 

In fact, Sun-like stars are the primary targets of 
space missions conducting planetary transit surveys like 
C0R0T and Kepler, where a large number of eclipsing bi- 
nary systems are monitored. We discuss the possibility 
of applying this method to bright eclipsing binaries mon- 
itored by Kepler, specifically the primary eclipse of low 
mass ratio binaries, where the low-amplitude signal can 
be detected by using data obtained during the lifetime 
of the mission. 
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We also show that there are in fact two photometric 
signals induced by the beaming effect during eclipse, one 
sensitive to the star's self rotation and spin-orbit angle, 
while the other is not. Both effects need to be accounted 
for in order to correctly interpret the light curve and 
estimate A. 

We describe the two effects in Section [5] and give ap- 
proximated analytic expressions for their amplitudes. In 
Section [3] we present a numerical simulation of an eclips- 
ing binary composed of an F and K type stars. We dis- 
cuss the results of the numerical simulation in Section 01 
including the possibility of detection with Kepler data, 
and possible degeneracies in the model. We summarize 
our work in Section [5] 

2. THE TWO BEAMING-INDUCED SIGNALS DURING 
ECLIPSE 

2.1. The photometric RM (PRM) effect 

The PRM effect results from the combination of two 
effects, the RM effect and the beaming effect, briefly de- 
scribed in the following paragraphs. 

During different phases of an eclipse, the eclipsing star 
blocks the light coming from different regions of the 
eclipsed star surface, that have different radial velocities 
(RVs) due to the rotation of the eclipsed star. The ob- 
served spectroscopic line profile is distorted and its center 
is shifted, resulting in the RM RV signal. The RM RV 
curve shape depends primarily on the sky-projected an- 
gle between the eclipsed star's spin and the system's or- 
bital angular momentum, A, the sky-projected rotational 
velocity of the eclipsed star, V rot sin/, and the secondary 
to primary radii ratio, r = R s /R v (e .g., lOhta et aT1l2005t 
IGimenezll2006t IGaudi fc Winnll2007l) . 

The beaming effect causes the observed flux of a light- 
emitting object to vary with its RV variation. When the 
object's RV varies periodically, as is the case in binary 
systems, it induces a periodic sinusoidal variation in the 
observed flu x. When the light of the secondary can be ne - 
glected (e.g. JMazeh fc Faiglerl2010HShporer eTaHl2(H(l . 
the amplitude of such a photometric modulation, in rel- 
ative flux, is otheam^K j 'c, where K is the RV amplitude 
of the primary and c the speed of light. The coefficient 
otbeam equals unity for bolometric light, and for a finite 
bandpass its value depends on the target's spectrum and 
the detector transmission curve. Since in stellar bina- 
ries the RV of both components modulates in opposite 
phase, the observed effect is the weighted difference of 
the individual bea ming effect amplitud es, weighted by 
their relative fluxes (jZucker et al.ll2007f) : 



^bcarn — 

C 



^ abeam ,p ** p * p 



abeam ,s ** s * s 



Fp + F s 



(1) 



where ab ea m,p, K p , and F p are the primary's beam- 
ing coefficient, RV amplitude and flux, respectively, 
while abeam, s, K s , and F s are the correspond- 
ing quantities for the secondary. Several authors 
have already observed t his effect, both from space 
e.g.. Ivan Kerkwiik et all I20T0I: IMazeh fc Faigien |20Tot 



iBloemenet al.l 1201 lHCarter etal.1 120111 ) and from the 
groun d (jShporer et al.l 120101 see also iMaxted et al.l 
12001 . 

When we apply the beaming effect to the eclipsed 
star rotational velocity, we get the photometric analog 



of the RM effect which we refer to as the photometric 
RM effect, or PRM. Photometrically, the blue shifted 
photons emitted from the star's hemisphere rotating to- 
wards the observer make that hemisphere appear slightly 
brighter than the other, redshifted hemisphere, rotating 
away from the observer. It is this apparent non uniform 
surface brightness that causes a photometric anomalous 
signal during eclipse. 

We adopt a coordinate system centered at the eclipsed 
star, and the Y axis directed towards the observer while 
the Y-Z plane contains the eclipsed star rotation axis. 
Therefore, the variation in relative flux during primary 
eclipse does not depend on the secondary position along 
the Y axis, but only on the position within the X-Z plane, 
and can be described by: 



A/ = Vilsm 



J J xl(x, z)dxdz 
J J T(x, z)dxdz 



(2) 



where $7 and T(x 1 z) are the eclipsed star rotation rate 
and surface brightness (not considering the PRM effect), 
respectively, and T> is the dilution factor, accounting for 
light from the eclipsing star, so A/ is the observed rel- 
ative flux difference. The integration includes the non 
eclipsed stellar surface and we assume flux of unity when 
the stars are out of eclipse. Eq. [5] is similar to the corre- 
sponding equatio n for the spectroscopic RM effect (e.g., 
lOhta et ail 120051 ). where the multiplication between fl 
and x gives the RV at position (x,z). 

An order of magnitude estimate of the PRM amplitude 
is: 



iPRM 



4a 



Vrot sin / 



beam " 



10" 



, V rot sin I 
10 km s- 1 0.1 



(3) 

We caution that Eq. [3] is only a rough estimate, and 
^4prm decreases for orbits with an inclination angle of 
i < 90 deg. The primary assumptions behind Eq. [3] are 
that the eclipsing star is relatively small and faint com- 
pared to the eclipsed star. 

2.2. The in-eclipse orbital beaming (InOrB) effect 

The orbital beaming effect follows the Keplerian mo- 
tion of the two stars. Therefore, at the time of con- 
junctions the orbital beaming modulation is an almost 
completely linear variation with time, especially for a 
circular orbit. When an eclipse occurs the observed frac- 
tional brightness of the two stars changes, leading to a 
variation in the weights in Eq. [TJ resulting in a deviation 
from the orbital beaming light curve. We refer to this 
deviation as the in-eclipse orbital beaming, or InOrB. 

The amplitude of this second effect depends on the 
variation of the eclipsed star's relative brightness, and its 
RV variation during eclipse. The latter can be approxi- 
mated in the following way, for a circular orbit. The RV 
amplitude of a star of mass M\ with a binary companion 
of mass Mi is: 



K x 



2na Mi sin i 
Po^b Mi + M 2 



(4) 



where a is the orbital semi major axis, P or b the orbital 
period, and i is the orbital inclination angle. The orbital 
phase at the time of eclipse start relative to the mid 
eclipse time is close to A^ w (i?i +i? 2 )Vl — b 2 /a, where 
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Ri and i?2 are the radii of the two stars, and b is the 
eclipse impact parameter^]. In addition, we assume the 
orbital phase covered during eclipse is a small part of the 
orbit, which results in: 



ARV C 



eclipse 



KiA()> 



2n(Ri +R 2 ) M 2 ahxi 
P~b Mi + M 2 



V1-& 2 



(5) 

where M\ and R\ are the mass and radius of the eclipsed 
star, respectively, and M% and R2 are the mass and radius 
of the eclipsing star, respectively. 

A rough estimate of the two effects amplitudes ratio 
gives: 



A 



InOrB 



-frr 



M 2 



Ri + i?2 sin? 



Aprm Porb M1 + M2 Rt sin/i ' 



(6) 



where P ro t,i is the eclipsed star rotation period. Eq. [6] 
shows that for certain systems the two effects can have 
similar amplitudes, especially for binaries of similar com- 
ponents where the stars' rotation is synchronized with 
the orbital period. 

As Fig. Q] and Fig. [2] show, for prograde orbits the In- 
OrB signal is opposite in sign to the PRM signal, and has 
a somewhat different shape, as the PRM signal at dif- 
ferent phases of the eclipse depends on the sky-projected 
distance of the eclipsing object from the eclipsed object's 
rotation axis, while the InOrB effect docs not. In fact, 
the InOrB effect does not depend on the eclipsed star 
rotation, and occurs even if the star did not rotate at all. 
As it does not hold information about spin-orbit align- 
ment, or any other new information, the InOrB effect is 
determined by parameters known from the eclipse light 
curve (considering only geometry and limb darkening) 
and the out-of-eclipse (i.e. orbital beaming) light curve. 

3. NUMERICAL SIMULATION 

To show the shape and amplitude of the light curves 
induced by the effects described in Section [2] we exam- 
ine a hypothetical eclipsing binary system and construct 
the light curves using numerical integration across the 
observed stellar sur faces. Our numerical co de is based 
on the code used in lShporer fc Brownl (|2011l . presenting 
a low-amplitude effect related to the spectroscopic RM 
effect). 

Our code assumes the two stars are rotating spheri- 
cal bodies, orbiting each other in a Keplerian orbit, and 
radiating as blackbodies with given effective tempera- 
tures. We further assume their surface brightness obeys a 
quadratic limb darkening law and we use limb darkening: 
coefficients based on the tabulation of lClaret fc Bloemenl 
(|2011ft , calculated specifically for the Kepler transmission 
curve. For simplicity we use a circular orbit. The beam- 
ing coefficient for each star was estimated by integrating 
the relevant blackbody spectrum across the Kepler trans- 
mission curve. 

The binary system we adopt here is composed of an 
F8V primary, with T e // jP =6250 K, and a KOV secondary, 
with T e // :S =5I50 K, orbiting each other every 10.0 days. 
We assumed solar metallically and an age of 1 Gyr, 
and der ived the stars' masses and radii using the Y 2 
models (|Yi et al.|[200l . We use M p = 1.20 M Q and 

5 defined as b = ^ C °V 



TABLE 1 

Parameters of the system used in the 

examples presented in flgs. 1 and 2. 



Parameter 


Value 


Units 


Orbital parameters: 






Orbital Period, P or t 


10.0 


day 


Eccentricity, e 





— 


primary stellar parameters: 






Mass, M p 


1.20 


M 


Radius, R p 


1.15 


-Rco 


Temperature, T e ff p 


6250 


K 


Rotational velocity, Vrot,p 


20 


km s 


primary stellar parameters: 






Mass, M a 


0.80 


Mq 


Radius, R s 


0.70 


Rr 


Temperature, T e ff a 


5150 


K 


Rotational velocity, V ro t,s 


5 


km s _1 



R1+R2 



Rp = 1.15 R Q for the primary, and M s = 0.80 M Q 
and R s — 0.70 Rq for the secondary. For the rotation 
rates we use V ro t, p — 20 km s _1 and V ro t,s — 5 km s _1 
for the prima r y and secondary, respectively, based on 
Mcibom ct al. (2011). All system parameters are listed 
in Table [TJ These parameters and the above assumptions 
imply an orbital semi major axis of 0.11 AU. 

Fig. [1] presents the various effects occurring during pri- 
mary (left panel) and secondary (right panel) eclipses, 
where we have assumed an edge-on orbit (inclination an- 
gle of i=90 deg, or an impact parameter of 6=0), and 
complete spin-orbit alignment for both stars (X p = A s = 
0). The eclipse duration for this system is 5.73 hour. 
The top row panels of Fig. Q] show the familiar eclipses 
light curves, in relative flux versus time, according to 
the geometry of the system and the limb darkening of 
the stellar surfaces. The effects induced by beaming 
are too small to be visually identified at the scale pre- 
sented in the top panels, so they are presented in the 
three bottom rows. The second to top row panels show 
the PRM light curve. For a prograde system, the ap- 
proaching side of the star is occulted first. Due to the 
beaming effect this side is slightly brighter than the re- 
ceding side, leading to a small flux decrease during the 
first half of the eclipse, and a flux increase during the sec- 
ond half. The schematic shape of the PRM light curve is 
similar to t he spectroscopic RM RV curve for stellar bi- 
naries ^(e^^^QbrecM^r^Ll [2003) and transiting planets 
(e.g., iGaudi fe Winn 2007). For the secondary eclipse, 
the PRM amplitude is much smaller than for the pri- 
mary eclipse, since the secondary is fainter, smaller, and 
rotates more slowly than the primary. Also, the PRM 
effect is seen only during ingress and egress of the sec- 
ondary eclipse (whose duration is 2.18 hours) since the 
star is completely occulted during the 1.37 hours of the 
full eclipse phase. The second to bottom row panels show 
the orbital beaming around the time of superior and in- 
ferior conjunctions. The linear trend seen in both panels 
is the orbital beaming, following the systems Keplerian 
motion. The small deviation from a straight line is the 
InOrB effect, described in Section l2~2l Finally, the bot- 
tom row panels show the two beaming-induced effects 
during eclipse, the PRM in red solid lines and the InOrB 
in blue solid lines where the linear orbital trend is sub- 
tracted, so these panels show only the anomalous anti- 
symmetric signals that result from the eclipsing nature 
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Fig. 1. — Photometric variability during primary (left) and secondary (right) eclipses of the F8V+K0V eclipsing binary whose parameters 
are listed in Table[TJ In this figure we assumed an edge-on orbit (i=90 deg) and spin-orbit alignment (A=0 deg) for both stars. From top 
to bottom: (a) eclipse light curve, accounting for geometry and limb darkening, in relative flux, (b) The PRM light curve, (c) The orbital 
beaming light curve, including the in-cclipse beaming, (d) The two beaming-induced effects in eclipse, PRM in red and the InOrB in blue. 
The dashed black line marks their sum. Table [TJ lists the parameters of the system analyzed here. As expected, during secondary eclipse 
the PRM effect is much smaller, since the secondary rotates more slowly and is fainter than the primary, and it occurs only during eclipse 
(occultation) ingress and egress. In fact, the PRM amplitude in the secondary eclipse is smaller than the InOrB amplitude. 
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Fig. 2. — Comparison of the beaming-induced signals during primary eclipse, for different impact parameters and spin-orbit angles. The 
left panel uses b=0 («=90 deg), and the right panel b=0.2 (i=89.14 deg). In both panels the red lines represent four different values of A, 
listed in the panels' legends. The blue solid line shows the InOrB light curve, which does not depend on A. 
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of the system. The dashed black line is the sum of the 
two effects, for reference. As can be seen in the bottom 
left panel, during primary eclipse the PRM amplitude is 
much larger than the InOrB amplitude. However, this 
relation is reversed in the secondary eclipse. 

In Fig. [2] we focus on the primary eclipse and show 
the beaming-induced signals for different eclipse impact 
parameters and spin-orbit orientations. The left panel 
shows the same system as in Fig. [T] and Table [T] but for 
four different values of Ai. The red lines show the PRM 
light curve for spin-orbit angles of (solid line), 45 (dot- 
dashed line), 90 (dotted line) and 180 deg (dashed line). 
The figure shows how the PRM signal decreases in am- 
plitude but remains antisymmetric for this completely 
edge-on case. The blue solid line shows the InOrB ef- 
fect, which is insensitive to spin-orbit angle, so it is the 
same for all four orientations. The right panel presents a 
similar system, only with an impact parameter of 6=0.2, 
which is equivalent to z=89.14 deg, and it shows how the 
PRM light curve becomes distorted for non-aligned and 
non-edge-on systems. 

4. DISCUSSION 

For the edge-on and aligned system presented in Fig.[TJ 
the PRM amplitude is 55 ppm for the primary and 2.2 
ppm for the secondary eclipses. These amplitudes are 
smaller than predicted by Eq. [3] due to dilution from 
the eclipsing star and also because that equation was 
derived by assuming the eclipsing object is much smaller 
than the eclipsed one. When this is not the case the 
integration across a significant fraction of the eclipsed 
stars' surface acts to decrease Aprm by averaging down 
the variation in the radial component of the rotational 
velocity. The InOrB effect amplitude for the primary is 
2.7 ppm, which is 5% of the PRM amplitude, not far 
from the 0.1 ratio predicted by Eq. [6] For the secondary 
eclipse the numerically measured InOrB amplitude is 6.2 
ppm, almost 3 times larger than the PRM amplitude, 
much larger than the ratio of 0.6 predicted by Eq. [51 
This is not surprising since the assumptions done when 
deriving Eq. [6] are invalid for the secondary eclipse. 

Both beaming-induced eclipse signals presented here 
are small, and can be detected only with ultra precise 
photometry. To estimate whether they can be detected 
by Kepler we assume an observed brightness of 12th mag- 
nitude, wher e according to the Kepler Mission sp ecifi- 
cations (e.g.. iBorucki et al.l l2010t Koch et "all 120101 ) the 
photometric precision is 20 ppm over 6.5 hours (to allow 
detection of the 13 hours long and 84 ppm deep transit 
of an Earth- like planet across a Sun- like star). Assuming 
white noise this translates to 72 ppm per 0.5 hour. Such 
a precision will of course not allow a detection from a 
single eclipse event even for the primary eclipse. Fortu- 
nately, Kepler's long lifetime, of 3.5 years, will allow to 
fold together many eclipses and decrease the scatter by 
y/N where N is the number of observed eclipse events. 
For the 10 day system adopted here it is fair to assume 
that 100 eclipse events will be observed, even when ac- 
counting for gaps following data download time or space- 
craft Safe Modes. Therefore, an accuracy of about 7 ppm 
per half an hour is predicted to be achieved over the life- 
time of the mission. Such a precision should be sufficient 
to detect the PRM effect during primary eclipse for a 
system similar to the one simulated here. 



Most of the Kepler targets are observed in the default 
Long Cadence mode. In this mode every 270 individ- 
ual consecutive exposures are stacked on-board, each of 
6.02 s exposure time and 0.52 s readout time. This re- 
sults in an effective cycle time of 29.4 min where photons 
are collected during 91.4% of that time. The high duty 
cycle allows us to treat each cycle as half an hour expo- 
sure. Hence, the individual eclipse events of the system 
shown in Fig.[T]are sampled with 12-13 exposures, which 
is sufficient for resolving the temporal flux variation al- 
though the variation across half an hour is averaged out. 
The smaller PRM signal during secondary eclipse makes 
it more difficult to detect, and probably not detectable 
with Kepler, although the amplitude increases for stars 
rotating more rapidly, which is expecte d for younger star s 
due to the Skumanich spin-down law (Skumanich 1972). 

The system we examined here has stellar components 
which are significantly different from each other. If we 
take them to be more similar, then the secondary eclipse 
signal becomes larger, although the primary smaller. For 
a system composed of two GOV stars with a self rota- 
tion velocity of 10 km s _1 , our numerical calculations 
show that the PRM amplitude will be close to 15 ppm. 
For such a system the InOrB amplitude will be 11 ppm, 
larger than for the case of the F8V+K0V system. 

The signal to noise ratio of the in-eclipse light curve 
is determined by the total amount of time the system 
spends in eclipse. Since for given stellar radii this scales 
with P -2 / 3 , eclipsing systems orbiting with a shorter or- 
bital period will require data taken over a shorter time 
span in order to reach the required precision. Taking the 
same 12th magnitude system adopted here, only with 
-Por6=l-5 day gives a precision of 7 ppm per half an 
hour after 150 days of data, within two Kepler quarters^. 
However, the eclipse duration of such a short period sys- 
tem is only 3.07 hours. This means that while Aprm 
remains the same, each eclipse will be covered by 6-7 
Long Cadence exposures, and each half of the eclipse by 
only 3-4 exposures. Moreover, when all Long Cadence 
measurements are folded together and binned the effec- 
tive averaging of the flux variability goes beyond half an 
hour. Therefore, using Long Cadence data to detect the 
PRM signal of short period eclipsing binaries will be dif- 
ficult, and even more so for distinguishing between the 
different possible spin-orbit orientations. Using Kepler 
Short Cadence data, with an effective integration time 
of 1 minute (9 readouts), could make it possible to de- 
tect the PRM effect in short period binaries. However, 
the Short Cadence mode is applied only to a small sample 
of 512 objects, about 0.3% of the Kepler targets, which 
are carefully chosen and observed in this mode usually 
during only one Kepler quarter. 

Measuring spin-orbit alignment for longer period sys- 
tems, with a period of 10 days, as simulated here, is also 
expected to be more useful to study binary formation, 
compared to short period systems. Processes resulting 
from tidal interaction between the two stars, including 
orbital circularization, synchronization between the or- 
bital p eriod and the stars' rotation, and spin-orbit align- 
ment (Mazeh 2008), have a time scale that is strongly 
dependent on the scaled semi- major axis, a/R p (e.g., 
iZahrt fl989h . The spin-orbit alignment specifically has 

6 Kepler quarters have a duration of about 90 days. 
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a relatively short time scale, two or three orders of mag- 
nitude shorter than the circularization time scale since 
it involve s a transfer of a smaller amount of angular mo- 
mentum (Mazeh 2008). Therefore, it is expected that 
for short period binaries the spin axes of the two stars 
had already aligned with the orbital momentum axis, and 
any initial missalignment — a relic from their formation 
history — had already faded away. Since for longer pe- 
riod systems, with a wider orbit, this process is slower, 
they are better targets for studying the initial alignment 
and the systems' formatio n process. For exam ple, the 
missaligned system DI Her (jAlbrecht et alJl2009t see also 
Mazeh 2008) has an orbital period of 10.6 day. 

We used the Kepler eclipsing binary catalog of 
iPrsa et all (|2011| ). including 1879 systems, to estimate 
the number of eclipsing binary systems whose photome- 
try will be accurate enough to allow looking for the PRM 
signal. We looked for systems whose eclipse duration is 
longer than 5 hours, and that given their orbital period 
and brightness (Kepler magnitude) the photometric S/N 
will be no less than half that of the 10 days and 12th 
mag system examined here. This includes 10 day sys- 
tems whose brightness is down to 13.5 mag, or 12 mag 
systems with period up to 40 days. This cut left about 
320 systems. When we restrict ourselves to low mass ra- 
tio systems, similar to the one we examined here, this 
number is about 50. A detection of the PRM effect and 
a measurement of the spin-orbit angle for the primary in 
even some of those systems will greatly extend the sam- 
ple of eclipsing binaries for which this angle is known. 

In our numerical analysis we ignored several low- 
amplitude effects, including ellipsoidal distortion due to 
tidal forces between the two stars and reflection of light, 
or heating of one star by the other. For a 10 day pe- 
riod binary both these effects are expected to have an 
amplitude along the orbit which is smaller th an that of 
the orbital beaming effect (|Zucker et al.| [2007). In addi- 
tion, during eclipses the ellipsoidal and reflection effects 
are at extrema, meaning their induced flux variation is 
minimal, while the orbital beaming flux is at its phase 
of fastest variation. The light curves of the ellipsoidal 
and reflection effects during eclipse are expected to be 
symmetric about the mid eclipse time, and, as for the 
orbital beaming, they can be predicted using modelin g 
of the out of eclipse light curve (|Faigler fc Mazeh~l l2011). 

We note that we ignore here the possibility of other 
sources of low level noise, such as stellar noise, meaning 
that this measurement requires quiet stars. In case the 
eclipsed star shows slight activity in the form of small 
spots on its surface, this can lead to sinusoidal-like flux 
variability at the rotational period. However, for the 
type of system simulated here the orbital and self ro- 
tation periods are not necessarily synchronized, so the 
spots-induced flux variability will be averaged out once 
folded on the orbital period, which is determined to high 
accuracy by the measurement of eclipses timing. 

Our method has several possible degeneracies, which 
if not carefully considered can lead to incorrect inter- 
pretation of the light curve, and a biased value for A. 
Significant eccentricity can produce asymmetric eclipse 
ligh t curves, depen ding on the argument of periastron, w 
(see iKippind 120081 for a description of the impact of ec- 
centricity on planetary transit light curves). Although, 
eccentricity can be estimated from the orbital beaming 



light curve, and, the value of ecosw is accurately deter- 
mined from the time difference between the primary and 
secondary eclipses. 

A second degeneracy results from the PRM signal be- 
ing mimicked by a small shift in mid eclipse time, which 
may give the false impression that the PRM signal does 
not exist in the light curve, suggesting a scenario where 
i = 90 deg and A = 90 deg, shown in Fig. [2] left panel. 
We looked into this possibility more closely by subtract- 
ing the eclipse light curve including a small shift in mid 
eclipse time, and without the PRM effect, from the light 
curve including the PRM signal and no time shift. The 
mid eclipse time shifts we examined were up to a few 
seconds, in sub second steps. We found that the peak 
to peak amplitudes of the difference light curves were 
at least 60 ppm, and the smallest difference is for time 
shifts close to 1 s. We expect Kepler photometry ob- 
tained during the time span discussed here to allow mid 
eclipse measurements of better than 1 s precision. There- 
fore, a solution which does not include the PRM effect (or 
assumes i — 90 deg and A = 90 deg) and as a result gives 
a wrong mid eclipse time due to the above mentioned 
degeneracy will have a significantly lower likelihood than 
a solution that does include the PRM effect. 

Modeling the PRM effect adds two parameters to the 
model, VrotsmI and A of the eclipsed star. An inde- 
pendent estimate of V TO t sin I can be achieved from a 
spectrum of the star, thereby constraining this param- 
eter so the information in the PRM signal is used for 
fitting only a single parameter (A). Moreover, an inde- 
pendent estimate of V ro t sin / is important for lifting a 
third degeneracy in our model, between A an d V rot sin 7, 
for edge-on orbits fe.g. JGaudi fe Winnll2007t ). 

Another possible degeneracy may result from the In- 
OrB signal being opposite in sign to the PRM signal (see 
Figures [T] and [2|), giving a false impression that the PRM 
signal is smaller in amplitude. In general, for non-aligned 
and non edge-on orientations, the combined light curve 
shape could be different than the PRM light curve, possi- 
bly leading to a biased estimate of A, even though for the 
primary eclipse shown in Figures [1] and [2] the PRM sig- 
nal is clearly dominant. However, as already noted, the 
InOrB light curve shape is insensitive to the spin-orbit 
angle or the star's rotation velocity. The only parameter 
to which both signals are sensitive to is the radii ratio 
r, which is accurately determined from the eclipse light 
curve even without accounting for the beaming-induced 
effects. Therefore, the InOrB light curve is determined 
independently of the PRM light curve, eliminating this 
possible degeneracy. 

Lastly, we note that the oblateness of stars, due to self 
rotation, may also affect eclipse light curves, especially 
when oblaten ess leads to significant gravity darkening. 
iBarnes! (|2009l ) showed how this will distort the shape of 
planetary transit light curves when the host star is a 
rapidly rotating early type star. However, since Sun-like 
stars are slow rotators this is expected to have a minor 
impact on the light curves (for the Sun for example the 
measured oblateness^ is 10~ 5 ). 

It is interesting to consider the application of the 
method presented here to a transiting planet orbiting 

7 The fractional difference in the radius at the equator and the 
poles. 
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a Sun-like star. However, as seen from Eq. [3J in this 
case Aprm is expected to be at the 1 ppm level, mak- 
ing it difficult, if at all possible to detect. The case of 
a transiti ng pla net orbiting an A type star is presented 
bv lGrootl (J2011I ). where the assumed rotational velocities 
are at the order of 100 km s _1 , and the presented light 
curve has an amplitude of about 10 ppm. 

5. SUMMARY 

We have presented the effects induced by beaming dur- 
ing eclipses of binary systems composed of Sun-like stars. 
Although a detection of these effects is clearly highly 
challenging, we have presented a numerical simulation of 
a F8V + KOV eclipsing binary where data accumulated 
over the Kepler Mission lifetime should allow to detect 
it during primary eclipse, and in turn measure the spin- 
orbit alignment of the primary. Given the known charac- 
teristic of Kepler eclipsing binaries we expect there will 
be about 50 systems similar to the system simulated here, 
where the PRM signal can be measured. In case the Ke- 
pler mission will be extended, this number will become 
significantly larger. So far the spin-orbit angle was mea- 
sured in only a few stellar binary systems, making its 
measurement in any additional system of high scientific 
value. 

The method presented here complements other meth- 
ods used for obliquity measurement in stellar binaries 
and star-planet systems. Those inc lude the measurement 
of the spe c troscopic RM eff e ct (1 Albrecht et all 120071 
[20091 l201ll : IWinnet al.l [20051 iTriaud et al.n20Tl . 



and 



the photometric methods presented by iBarned (|2009|) . 
for rapidly rotating e arly t ype planet host stars, and by 
ISanchis-Oieda et al.l (120111 ). for active planet host stars. 
In addition. iGrootl (|2011| ) described the application of 
this method to early type eclipsing stellar binaries, a 
transiting star-planet system with an early type host, 
and eclipsing double white dwarfs. 

Pursuing the measurement of the spin-orbit angle in 
eclipsing binaries with Sun-like components using pho- 
tometry, as we described here, has a few advantages rel- 
ative to other kind of systems and other methods. First, 
Kepler is already collecting high precision photometry 
and is targeting primarily Sun-like stars. Kepler's data 
will eventually have a long time span, and will become 
publicly available. Second, the deep spectral absorption 
lines of late type stars allow obtaining high precision RVs, 
where the photometric results can, in principle, be con- 
firmed. A similar confirmation will be more difficult for 
early type stars, due to the nature of their spectra, and 
even more so for double white dwarfs, since the eclipses 
are too short to allow sufficient exposure time. 
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